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ABSTRACT 

We use observations of three galaxies hosting low-ionization nuclear emission-line regions (LINERs; 
NGC 404, NGC 4736, and NGC 4579) with the Chandra X-ray Observatory to study their power sources. 
We find very diverse properties within this small group: NGC 404 has an X-ray-faint nucleus with a soft, 
thermal spectrum, NGC 4736 harbors a plethora of discrete X-ray sources in and around its nucleus, 
and NGC 4579 has a dominant nuclear point source embedded in a very extended, diffuse nebulosity. 
From their multi-wavelength properties we conclude the following about the power sources in these 
LINERs: the nucleus of NGC 404 is the site of a weak, compact starburst, whose X-ray emission is due 
Qh to gas heated by stellar winds and supernovae, NGC 4736 appears to be in a recent or aging starburst 

phase, where the X-ray emission is dominated by a dense cluster of X-ray binaries, and NGC 4579 
is powered by accretion onto a supermassive black hole. We detect 39 discrete sources in NGC 4736 
and 21 in NGC 4579, most with L x > 10 37 erg s" 1 . One source in the disk of NGC 4579 appears 
to be an ultraluminous X-ray binary with Lx(2 — 10 keV) = 9 x 10 39 erg s _1 , but it could also be a 
background quasar. The most luminous discrete sources have simple power-law spectra, which along 
with their luminosities suggest that these are X-ray binaries accreting near or above the Eddington 
rate for a neutron star. By comparing the luminosity functions of discrete X-ray sources in these and 
other galaxies we find a potential connection between the age of the stellar population and the slope of 
the cumulative X-ray source luminosity function: galaxies with primarily old stellar populations have 
steeper luminosity functions than starburst galaxies. We suggest that this difference results from the 

■ contribution of high-mass X-ray binaries from the young stellar population to the upper end of the 
luminosity function. 
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far-reaching implications. The frequency of LINERs that 
are weak AGNs sets a lower limit to the frequency of su- 
permassive black holes in the centers of galaxies. If the 
majority are really accretion-powered sources, they would 
be the most common type of AGN in the present-day 
universe, representing the faint end of the AGN luminos- 
ity function. They could be the remnants of once fero- 
cious quasars, and as such their demographics constrain 
the endpoint of quasar evolution. Moreover, the observed 
spectral energy distributions suggest that the accretion 
flows of weak AGNs in LINERs have a different struc- 
ture than in Seyferts (e.g., Lasota et al. 1996; Ho 1999). 
They may be described by advection dominated accre- 
tion flows (ADAFs; Narayan & Yi 1994, 1995), adiabatic 
inflow-outflow solutions (ADIOS; Blandford & Begelman 

1999) , or related systems, which are interesting in their 
own right for understanding the growth of supermassive 
black holes and their radiative efficiency In those LINERs 
powered by compact starbursts, the very presence of the 
associated young star clusters, and also their properties, 
constrain formation scenarios of galaxy bulges (e.g., Bokcr 

2000) . More specifically, the tidal influence of a compact 
star cluster in a galaxy center can destroy a pre-existing 
galactic bar that led to the formation of the bulge in the 
first place (Carollo 1999). 

In view of the importance of LINERs as a class, we un- 
dertook a study of their X-ray properties with the Chandra 
X-Ray Observatory, in order to address some of the above 
issues. The results of previous X-ray observations have 
left ambiguities since older instruments had either high 
spatial resolution (e.g., the Einstein and ROSAT HRIs) 
or spectroscopic capabilities (e.g., the Einstein IPC, the 
ROSAT PSPC, and the ASCA SIS and GIS) but not both. 
The combination of high spatial resolution and hard X- 
ray response afforded by Chandra makes it possible to re- 
solve such ambiguities by putting nearby LINERs under 
the microscope and examining their power sources in great 
detail. Accordingly, we have adopted the strategy of ob- 
taining long exposures of a small number of LINERs. Our 
program thus complements shallower Chandra surveys of 
large, statistically complete samples of nearby galaxies, 
which can address different types of questions (e.g., Ho et 
al. 2001). 

In this paper we present the results of our detailed study 
of 3 LINERs observed during Chandra Cycle 1. Included 
are two objects observed as part of our own guest-observer 
program (NGC 4736 and NGC 4579) and one object ob- 
served as part of a different program by a different group 
(NGC 404; the data were obtained from the Chandra pub- 
lic archive; program 01700766, P.I. P. Lira). We present 
the properties of the target galaxies in §2 and give the de- 
tails of the observations and data reduction in §3. After 
describing the X-ray morphologies of the three galaxies in 
§4, we embark on a detailed analysis of their X-ray sources 
in §5-7. We conclude with a discussion of the results and 
their implications for the nature of LINERs as a class in 
§8 and §9. 

2. TARGETS AND THEIR PROPERTIES 

Our targets are nearby, well-studied galaxies. They are 
listed in Table 1 in order of increasing distance, along 
with their basic properties. NGC 404 is a small elliptical 



galaxy while NGC 4736 and NGC 4579 are nearly face-on 
spirals (the latter is a barred spiral). All three galaxies 
were observed in the Hubble Space Telescope (HST) ul- 
traviolet snapshot survey of nearby galaxies (Maoz et al. 
1995, 1996) and were found to harbor unresolved, nuclear 
point sources, hinting to the presence of an AGN. Later 
UV spectroscopy of NGC 4579 and NGC 404 (also with 
the HST) showed the former to be a bona fide AGN with 
broad emission lines. The UV spectrum of NGC 404, on 
the other hand, resembles very closely that of a starburst, 
with prominent absorption lines from hot stars and no in- 
dication of an AGN. This result also confounds the case 
for an AGN in NGC 4736, since it shows that a point- 
like UV source could well be associated with hot stars. 
Narrow-band Ha and [O III] images of the three galax- 
ies, obtained from the ground and with the HST (Pogge 
1989; Pogge et al. 2000), resolve their nuclear emission- 
line regions down to scales of less than 10 pc and show 
complex filamentary structure. It is noteworthy that the 
Ha-emitting region of NGC 4579 is easily resolved from 
the ground; it extends along the direction of the bar of the 
host galaxy, with a projected size of approximately 3 kpc. 
The nucleus of NGC 4736 is surrounded by a star-forming 
ring of diameter 2 kpc, seen very clearly in the Ha images 
of Pogge (1989), as well as in the radio continuum image 
of Turner & Ho (1994). 

The target galaxies have been observed several times in 
the X-ray band. The two spirals are fairly bright X-ray 
sources in both soft and hard X-rays. Their ASCA spec- 
tra (Terashima et al. 1998; Roberts, Warwick, & Ohashi 
1999) can be described by a composite model consisting 
of a hard power law with a Seyfert-like photon index of 
~ 1.7 plus thermal plasma emission with a characteristic 
temperature in the range 0.1-1 keV, just like the spec- 
tra of most LINERs (Ptak et al. 1999). The observed 
luminosities are Lx(2 — 10 keV) = 1.5 x 10 41 erg s^ 1 
for NGC 4579 and L x (2 - 10 keV) = 3 x 10 39 erg s" 1 
for NGC 4736. Moreover, the ROSAT HRI image of 
NGC 4736 shows remarkable extended structure in the 
form of a halo surrounding the nucleus (Halderson et al. 
2001), which is comparable in size with the star-forming 
ring seen in the optical emission-line and radio continuum 
images. The nucleus of NGC 4736 was even resolved by 
the ROSAT PSPC, which measured a soft X-ray lumi- 
nosity of L x (0.1 - 2 keV) = 3 x 10 39 erg s^ 1 (Cui et al. 
1997). In marked contrast, NGC 404 is extremely faint 
even though it is the nearest of the three targets. It was 
barely detected in soft X-rays by the ROSAT HRI with a 
luminosity of L x (0.1 - 2.4 keV) w 5 x 10 37 erg s" 1 (Ko- 
mossa et al. 1999), while an ASCA observation yielded 
only an upper limit on the hard X-ray luminosity of 
L x (2-10 keV) < 5x 10 37 erg s" 1 (Terashima et al. 2000a). 

3. OBSERVATIONS AND DATA SCREENING 

The three target galaxies were observed with the Ad- 
vanced CCD Imaging Spectrometer (ACIS; Garmire et al. 
2001, in preparation) on Chandra in 1999 December and 
2000 May. The instrument was operated in faint mode. 
A log of observations, which includes the observation date 
and exposure time for each target, is given in Table 2. 
Their nuclei were centered on the aimpoint of the S3 CCD, 
which had not suffered significant radiation damage dur- 
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ing the early stages of the mission. In the observations of 
NGC 4736 and NGC 4579 only one quarter of the S3 CCD 
was read out; all other CCDs were turned off in order to 
reduce the frame time to 0.84 s, and thus to reduce the 
probability of photon pileup in the central pixels of bright 
point sources. Even with this measure, however, the bright 
point source in the nucleus of NGC 4579 still suffered from 
pileup, which complicates the spectral analysis (the cir- 
cumnuclear regions of NGC 4579 were not affected). In 
the observation of NGC 404, 6 CCDs were turned on (4 
from the ACIS-S array and 2 from the ACIS-I array), re- 
sulting in a frame time of 3.2 s. We consider here only 
the part of the image included on the S3 chip since this 
encompasses the entire optical galaxy. 

The initial data screening was carried out with the CIA0 
software package, developed and distributed by the Chan- 
dra X-Ray Center (CXC). Since this software was undergo- 
ing revisions over the course of our analysis of the data, we 
processed the data using both versions 1.0 and 2.0 to check 
for consistency. Some of the later steps in the data analy- 
sis were also checked with version 2.1 of CIA0. In summary, 
the screening consisted of selection of events with grades 
0, 2, 3, 4, and 6 and exclusion of events recorded dur- 
ing times of poor aspect solution and background flares. 
In the process we also removed the 0'.'5 pixel randomiza- 
tion introduced during the initial data processing, which 
amounts to an unnecessary smearing of the image. After 
screening, we produced images and spectra from the pho- 
ton event lists using the CIA0 tools and the TARA software 
package (Broos et al. 2000). Below we discuss the mor- 
phology of the target galaxies as seen in the X-ray images 
and compare with what is seen in other bands. 

4. MORPHOLOGY 

Figure 2 shows a montage of the 0.5-10 keV Chandra 
images of the target galaxies. For each galaxy we present 
an image encompassing most of the ACIS-S3 field of view 
as well as a zoomed display of the nucleus. These im- 
ages show exquisite detail, which we describe below in a 
separate paragraph for each galaxy. In our comparison 
with images taken at other wavelengths we have assumed 
that the "nucleus" of each galaxy has the same coordi- 
nates at all wavelengths. This is certainly true within the 
a priori relative astrometric uncertainties between images, 
which are approximately 1"5 (dominated by the Chandra 
aspect solution; Hornschemeier et al. 2001). Because of 
the small field of view and bandpasses of the Chandra and 
HST images that we are comparing, we are not able to 
find matching sources that can serve to improve the rela- 
tive registration of the images. We have also searched the 
USNO-A2.0 catalog (Monet et al. 1996) for optical coun- 
terparts to the detected X-ray sources, but found none. 

• NGC 404- ~ The only significant X-ray emission 
within the optical extent of NGC 404 comes from 
a compact region associated with the nucleus. Fig- 
ure 2 (top right panel) shows an image of the nu- 
cleus deconvolved via the Lucy-Richardson algo- 
rithm. There is clearly-resolved extended emission 
to the west and southwest of the nucleus, which af- 
ter deconvolution takes the form of discrete "blobs." 
The extended X-ray emission coincides with the ex- 



tended UV and line emission seen in the HST images 
of NGC 404, although there is also some low-level 
X-ray emission as far out as 10" from the center. 
Most of the flux is concentrated within 1" from the 
nucleus, which corresponds to a linear scale of 12 pc. 
The X-ray spectrum is extremely soft, as we show 
in later sections. This fact along with the shell- 
like structure of the X-ray and emission-line images 
suggests that the X-ray source is hot gas, very likely 
associated with a superbubble blown out by a com- 
pact starburst in the nucleus of NGC 404. Contrary 
to the two spiral galaxies in our collection and other 
galaxies observed with Chandra (including both spi- 
rals and ellipticals; Sarazin, Irwin, & Bregman 2001; 
Blanton, Sarazin, & Irwin 2001; Ward et al. 2000; 
Mushotzky et al. 2000; Angclini, Lowenstein, & 
Mushotzky 2001; Lira et al. 2001; Fabbiano, Zezas, 
& Murray 2001) we do not find any discrete X-ray 
sources away from the nucleus of NGC 404, down to 
a limiting luminosity of 1 x 10 36 erg s _1 . (Although 
there are many other sources in the field of view, 
none of them fall within the effective radius of the 
galaxy reported by Faber et al. 1989). 

• NGC 4736. - The image of NGC 4736 shows a 
plethora of discrete X-ray sources, with a concen- 
tration around the nucleus of the galaxy. The col- 
lective emission from these sources, which are most 
likely X-ray binaries, makes up the extended, unre- 
solved "halo" seen in the RO SAT SRI image, and 
accounts for most of the X-ray flux from this galaxy. 
In the innermost 400 x 400 pc of the galaxy we find a 
very dense cluster of 10 discrete sources embedded 
in a diffuse but clumpy halo. The four brightest 
sources (3 of which are in the central cluster) are la- 
beled as X-l through X-4 on the image (in order of 
decreasing X-ray brightness). We find that within 
the combined astrometric uncertainties of Chandra 
and the HST, the southern of the two UV point 
sources in the image of Maoz et al. (1995) coincides 
with source X-2, while the northern of the two UV 
sources appears to have no X-ray counterpart. The 
detection of large numbers of discrete X-ray sources 
both in NGC 4736 and in NGC 4579 (below) is un- 
doubtedly aided by the fact that the two galaxies are 
oriented face-on and have a low Galactic absorption 
column. The spectra of the discrete X-ray sources 
are fairly hard and their luminosities reach values 
well over 10 38 erg s _1 . The spectrum of the diffuse 
emission, on the other hand, is softer with a total lu- 
minosity comparable to that of one of the brightest 
discrete sources. We describe the spectral proper- 
ties of the discrete sources and the diffuse emission 
in more detail in later sections. The large concen- 
tration of discrete X-ray sources may be associated 
with a recent episode of star formation, which is still 
ongoing to some degree in the circumnuclear star- 
forming ring. We discuss this possibility further in 
§8.2. 

• NGC 4 579. - The X-ray emission from NGC 4579 is 
dominated by an unresolved source coincident with 
the nucleus, which contributes 90% of the 2-10 keV 
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X-ray flux. Wc identify this source with the AGN, 
which is the source of the broad optical and UV 
emission lines (Barth et al. 2001; Maoz et al. 1998). 
The AGN is embedded in an extended, diffuse X- 
ray halo, which is elliptical in shape. The halo ma- 
jor axis is ~ 2 kpc in length, aligned approximately 
parallel to the bar of the host galaxy. The X-ray 
halo appears to coincide exactly with the Ha emis- 
sion region seen in the narrow-band image of Poggc 
(1989), while the much more compact Ha-emitting 
region seen in the HST images of Pogge et al. (2000) 
corresponds to the innermost, brightest part of the 
X-ray halo, immediately surrounding the AGN. The 
image of NGC 4579 also shows about 20 discrete X- 
ray sources scattered throughout the disk of the host 
galaxy. These are presumably X-ray binaries, just 
as in NGC 4736. One of these sources, located at 
a projected distance of 8.7 kpc from the AGN at a 
position angle of 206°, is an ultraluminous X-ray bi- 
nary (UXB), with L x (2-10 kcV) = 9xl0 39 erg s" 1 . 
It is unresolved by Chandra, which sets an upper 
limit of 120 pc to its diameter. 

In the next section we present fits to the spectra of bright 
sources found in the three galaxies, while in later sections 
we study their variability and population properties. 

5. MODEL FITS TO THE SPECTRA OF BRIGHT SOURCES 
5.1. Extraction of Spectra 

To study the spectral properties of X-ray sources in the 
target galaxies, we extracted their spectra using CIA0. 
In particular, we extracted spectra of (a) the nucleus of 
NGC 404 and patches of diffuse emission in the nuclei of 
NGC 4736 and NGC 4579, (b) the AGN in NGC 4579, 
and (c) the 4 brightest discrete sources (X-l through X- 
4) in NGC 4736 and the cxtranuclear UXB (i.e. X-l) in 
NGC 4579. In all cases the spectra were extracted from 
a circular region centered on the source. The radius of 
this region was generally 4.1 pixels (corresponding to 2") 
for discrete sources, which encompasses more than 95% 
of the flux from a point source. In the case of the AGN 
in NGC 4579, which is quite bright and has a relatively 
hard spectrum, we enlarged the radius of the extraction 
region to 6.2 pixels (or 3") to include the wings of the 
point-spread function (PSF). In the case of the extended 
emission in NGC 4736, which underlies a crowded field 
of discrete sources, we extracted a spectrum from a re- 
gion that falls SE of source X-2 and SW of source X-l. 
This region has a radius equal to a discrete source extrac- 
tion radius and is free of discrete sources 5 . In the case 
of the extended emission from NGC 4579, we extracted 
a spectrum from an annulus centered on the AGN, with 
an outer radius of 20 pixels (or 9"8). This annulus ex- 
cludes the AGN and has holes to exclude discrete sources 
that fall within it, as well as the "readout streak". Its 
outer radius is small enough not to cross to a region of the 
CCD read out by a different amplifier. The response ma- 
trix and effective area curves corresponding to each source 
were computed automatically by CIA0 after the spectrum 
was extracted. 

A spectrum of the background for each observation was 

5 The absence of discrete sources was verified both by visual inspects 



extracted from a region of radius 50 pixels (or 24'.'5) cen- 
tered in a part of the field of view that appears free of 
discrete sources or diffuse emission. The 0.5-10 keV back- 
ground count rates in NGC 404, NGC 4736, and NGC 4579 
are, respectively, 6.4 x 10~ 7 ,2.5 x 10~ 6 , and 1.3 x 
10~ 6 s _1 pixel -1 , implying a total of 0.9, 6.1, and 2.3 
background counts per point source, respectively (assum- 
ing a radius of 2"). This background determination, how- 
ever, applies to the particle and diffuse sky background 
only. The effective background for the discrete sources of 
interest can be considerably higher (and with a different 
spectrum) because these sources are embedded in diffuse 
emission from the host galaxy. A serious, additional com- 
plication arises from the fact that the diffuse emission is 
patchy and uneven, which prohibits us from using its spec- 
trum, normalized by a geometric factor, as an effective 
background spectrum. To overcome this difficulty, wc in- 
clude the effective background in the model that we fit to 
discrete source spectra. More specifically, we fit the spec- 
trum of the diffuse emission first and the model we derive 
for it is added to the spectral model for each source. All 
the parameters of the effective background model are held 
fixed, except for a normalization factor. 

The spectra of the diffuse sources in NGC 4736 and 
NGC 4579 could also be contaminated by photons from 
the wings of the PSF of nearby, bright discrete sources. 
These contaminating photons have a relatively hard spec- 
trum since the PSF is more extended at higher energies. 
To eliminate this contamination we simulated the images 
of point sources that could affect the spectra of the dif- 
fuse sources in the two galaxies using the MARX tool. From 
the simulated images we extracted the spectra of contam- 
inating photons and subtracted them from the observed 
spectra of diffuse sources. 

Spectral fits were carried out using XSPEC v. 11. 0.1 (Ar- 
naud 1996). The spectra were truncated at the low end at 
0.5 keV and at the high end at energies generally between 
2 and 7 keV, depending on the signal-to- noise ratio (S/N). 
The spectral models described below were modified by in- 
terstellar photoelectric absorption with cross-sections as 
given by Morrison & McCammon (1983). The absorbing 
columns were constrained to be at least as large as the 
Galactic columns given in Table 1. The abundance ratios 
used in thermal plasma emission models are those of An- 
ders & Grevesse (1989). The error bars given on model 
parameters correspond to 90% confidence intervals (for 2 
interesting parameters), unless specified otherwise. 



5.2. Spectra of Diffuse Sources 

To model the spectrum of the diffuse emission in 
NGC 4736 and NGC 4579, as well as the spectrum of the 
nucleus of NGC 404, we tried combinations involving ther- 
mal plasma emission components (Mewe, Gronenschild, 
& van den Oord 1985; Mewe, Lemen, & van den Oord 
1986; Kaastra 1992; Liedahl, Osterheld, & Goldstein 1995) 
with and without a power law. Single-temperature plasma 
models provide acceptable fits to the spectra of NGC 404 
and NGC 4736, while a 3-temperature model is required 
for NGC 4579. The parameters of the best-fitting models 
are summarized in Table 3, and the spectra themselves are 

and by the systematic search described in §7 
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shown in Figure 3 with models superposed. The tempera- 
tures for the first two objects are around 0.6 keV, similar 
to what was inferred from previous ASCA observations of 
many LINERs (e.g., Ptak et al. 1999). In NGC 4579 one 
of the thermal emission components also has a tempera- 
ture of 0.6 keV, while the other two temperatures bracket 
this value. It is likely that there is a range of temper- 
atures in all cases, although this is only discernible in 
NGC 4579 where the S/N is highest. It is also notewor- 
thy that the heavy-element abundances are considerably 
below their Solar values in NGC 404. In NGC 4736 and 
NGC 4579 the abundances are consistent with the Solar 
values. Keeping the abundance value in the latter object 
fixed to Solar yields slightly different temperatures than if 
the abundances were free to vary, as we indicate in Table 3. 

Models consisting of thermal plasma emission plus a 
power-law can also yield acceptable fits, although such 
models are not preferable to the ones described above in 
a strict statistical sense. In the case of NGC 404, the 
power-law component has a photon index of 2.5+J'g, a 0.5- 
2 keV luminosity of 7 x 10 36 erg s _1 and an extrapolated 
2-10 keV luminosity of 5 x 10 36 erg s" 1 , equivalent to that 
of a single X-ray binary. In NGC 4579, the hottest thermal 
component can be substituted with a power-law compo- 
nent of photon index 1.5*qj5 (f° r % = Z Q ). More impor- 
tantly, however, the addition of a power-law to the model 
for NGC 404 allows us to fit the spectrum adequately by 
keeping the abundances fixed to their Solar values. The 
fundamental reason for this is that the emission-line com- 
plexes in the observed spectrum appear to be weaker than 
what thermal plasma models predict. Thus the addition 
of the featureless power law dilutes the line equivalent 
widths, and allows for acceptable fits with Solar metal 
abundances. With pure thermal plasma models, on the 
other hand, we can only achieve acceptable fits with metal 
abundances significantly less than Solar. 

There are two ways of interpreting the dilution of the 
thermal plasma lines of NGC 404. One possibility is 
that the abundances are indeed lower that Solar, while 
the other possibility is that there is a hard X-ray source 
with a featureless power-law spectrum embedded in its nu- 
cleus. As we discuss in more detail in §8.1, the nucleus of 
NGC 404 is best interpreted as a mild starburst, which 
most likely hosts newly formed X-ray binaries. Therefore, 
we favor the latter explanation of the apparent weakness 
of the thermal plasma emission lines. 

5.3. Spectra of Discrete Sources 

We fitted the spectra of bright discrete sources with 
a variety of models inspired by the results of studies 
of X-ray binaries in the Milky Way and the Magellanic 
Clouds. The X-ray spectra of Galactic X-ray binaries 
are typically described by composite models consisting 
of cither a multi-temperature accretion disk spectrum 
plus a Comptonizcd blackbody spectrum (motivated by 
TENMA observations; Mitsuda et al. 1984, 1989) or by 
a single-temperature blackbody spectrum plus unsatu- 
rated Comptonization (motivated by EXOSAT observa- 
tions; White, Stella, & Parmar 1988). More recent ob- 
servations with BeppoSAX and RXTE confirm these find- 
ings and show that these models can sometimes be indis- 
tinguishable, even with high- 5/ AT spectra (see review by 



Barret 2001, and references therein). The above models 
apply to objects in their "high" luminosity state (Lx ~ 
10 37 — 10 38 erg s _1 ). As the luminosity approaches and 
exceeds 10 38 erg s _1 , the spectral shape resembles a sim- 
ple power law very closely. Accordingly, we experimented 
with combinations of soft and hard spectral components, 
with the soft component taken to be either a black body 
or a multi-color disk and the hard component taken to be 
a Comptonized blackbody spectrum (Titarchuk 1994) or 
a simple power law. We also experimented with single- 
component models made up of either a simple power law 
or a Comptonized blackbody spectrum. 

We find that all spectra can be adequately described by 
a simple power law (resulting in xt < 1-1) , with photon 
indices in the range 1.1-1.8 and 2-10 keV luminosities be- 
tween 4 x 10 38 and 9 x 10 39 erg s _1 (the most luminous 
source is, of course, the UXB in NGC 4579). The param- 
eters of the individual power-law models are given in Ta- 
ble 4, while the observed spectra with models superposed 
are shown in Figure 4. In two cases we measure absorbing 
columns that are almost an order of magnitude higher than 
the Galactic column, suggesting significant absorption in 
the vicinity of the source. These large columns should be 
regarded with caution, however, because they do not ap- 
pear to be so large in all the models that we experimented 
with. The fact that simple power-law models provide a 
good description of the observed spectra is consistent with 
the very high luminosities of these objects, just as in Galac- 
tic X-ray binaries. Most of the other spectral models de- 
scribed above fit the observed spectra equally well. How- 
ever, their larger number of free parameters makes them 
less appealing in a strict, statistical sense. Moreover, de- 
generacy between the model parameters and the limited 
S/N and spectral coverage prevent us from placing mean- 
ingful constraints on the parameter values. We note, nev- 
ertheless, that the fits of simple Comptonization models 
yield parameters comparable to those found in Galactic 
X-ray binaries, namely, seed blackbody temperatures be- 
tween 0.2 and 1 keV, Comptonizing plasma temperatures 
between 2 and 5 keV, and optical depths between 2 and 
6. The UXB in NGC 4579 is a notable exception in that 
a two-component model consisting of a multi-color disk 
(fc7i nnor = 0.4 ± 0.2 keV) and Comptonized blackbody 
(kT hh = 0.5±g;| keV, kT compt < 200 keV, r < 1) provides 
a somewhat better fit than a simple power-law model (% 2 
decreases by 4.71 at the expense of adding 3 more free 
parameters, which is significant only at the 73% level). 

5.4. The spectrum of the AGN in NGC 4579 

Because of its high (observed) count rate of 0.719 ± 
0.005 s-\ the spectrum of the AGN in NGC 4579 is 
severely distorted by the effects of pileup. In practice, 
the main consequence of pileup is to make the spectrum 
appear harder than it really is. This is a result of two ef- 
fects: (a) two or more lower-energy photons arriving in the 
same pixel within one frame time are misidentified as sin- 
gle, higher-energy photons, and (b) two or more photons 
arriving in the same or adjacent pixels in one frame time 
are completely missed because they produce non-standard 
event grades (this affects soft photons more than hard pho- 
tons because of their greater numbers). Pileup cannot 
yet be incorporated into the detector response matrices 
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and effective area curves, thus piled-up spectra must be 
modeled using special techniques. To this end, we have 
used the simulator-based, forward-fitting tool LYNX, devel- 
oped by the ACIS instrument team (Chartas et al. 2000). 
This tool simulates the distribution of counts per energy 
channel for an assumed input spectrum and iteratively 
refines the parameters of the model until the simulated 
count distribution matches the observed one. The sim- 
ulation takes into account photon trajectories (with the 
help of the MARX tool) and pileup. Uncertainties in the 
model parameters are computed statistically from the dis- 
tribution of their values over a large number of realiza- 
tions of the model. Using the above method, we were 
able to fit the AGN with a simple power-law model, with 
the absorbing column density fixed at the Galactic value. 
The resulting photon index is T = 1.88 ± 0.03 and the 
2-10 keV flux and luminosity (with pileup taken into ac- 
count) are, respectively, Fx = 5.2 x 10~ 12 erg cm~ 2 s _1 
and Lx = 1.7 x 10 41 erg s _1 . This model yields a reduced 
X 2 value of 0.773 for 646 degrees of freedom. The observed 
spectrum is shown in Figure 5 with the best-fitting model 
superposed for comparison. We note for reference that fit- 
ting the observed spectrum without taking pileup into ac- 
count yields T = 1.20 and Fx — 6.3 x 10 -12 erg cm~ 2 s _1 . 
The degree of pileup can be quantified by the "pileup frac- 
tion" , which is defined as the ratio of the excess counts in a 
selected energy window to the true counts in that window. 
In this case we estimate the pileup fraction in the 0.5- 
10 keV band to be / p (0.5- 10 keV) = -0.51 and the pileup 
fraction in the 6-7 keV window to be / p (6- 7 keV) = 0.33. 
In other words, the 0.5-10 keV count rate is 1.51 times less 
than the true value and the and the 6-7 keV count rate is 
1.33 times more than the true value. 

To check the results of the pileup simulation, we ex- 
tracted spectra of the AGN from an annulus that includes 
only the wings of the PSF and excludes the core where the 
pileup occurs. We took the outer radius of the annulus to 
be 6" and we experimented with inner radii in the range 
(f.'5 — l'/2. To fit these spectra we modified the effective 
area curves by an energy-dependent scale factor that takes 
into account the fact that the size of the PSF depends on 
energy (this factor was determined with the help of syn- 
thetic PSFs simulated with the MARX tool). We found that 
these spectra became progressively steeper as the inner ra- 
dius of the annulus was increased, i.e., as more of the core 
of the PSF was excluded. For inner radii greater than 1" 
the photon index settled to a value consistent with the 
LYNX result. 

We also searched for the Fc Ka line detected in previous 
ASCA observations in 1995 and 1998 by Terashima et al. 
(1998, 2000c), but we did not detect it. The method we 
adopted was to isolate the spectrum in the 4-8 keV range 
and fit it with a model consisting of a power-law contin- 
uum plus Gaussian emission lines at 6.4, 6.7, and 6.9 keV. 
The power-law model for the continuum was only a pa- 
rameterization, with no particular physical significance. 
The effect of pileup on this portion of the spectrum is to 
raise the continuum level and reduce the photon index by 
transforming pairs of lower-energy photons into apparent 
high-energy photons. Thus, it effectively "dilutes" the ob- 
served equivalent width of any lines in this window. The 
amount of dilution of lines between 6 and 7 keV is the 



pileup fraction estimated above, / p (6 — 7 keV) = 0.33. 
Thus all equivalent widths or upper limits determined by 
our procedure can be rectified by multiplying by 1 + / p . 
Accordingly, we plot in Figure 6 the 99%-confidence up- 
per limits to the line equivalent width as a function of the 
assumed velocity width (FWHM). We also mark the mea- 
sured equivalent widths and FWHM from the older, ASCA 
observations, which appear to be inconsistent with our up- 
per limits. It is noteworthy that the photon index and 
2-10 keV flux that we measure, also differ somewhat from 
what was measured by ASCA in 1995 (r = 1.72 ± 0.02, 
Fx = 4.3 x 10 -12 erg cm~ 2 s _1 ), but are consistent with 
the ASCA measurements from 1998 (r = 1.81 ± 0.06, 
Fk = 5.3 x 10~ 12 erg cm~ 2 s _1 ). Although it is conceiv- 
able that the Fe Ka line intensity has varied between the 
the Chandra and ASCA observations, we regard this as a 
rather unlikely explanation since the luminosity and spec- 
tral index appear to be the same now as in the 1998 ASCA 
observation. 

A more plausible and likely explanation for the differ- 
ence between the ASCA and Chandra measurements is the 
fact that the ASCA spectra were extracted from a much 
larger region than the Chandra spectra and include emis- 
sion from diffuse thermal plasma as well as many discrete 
sources. To put this difference in perspective, we note that 
for the AGN in NGC 4579 the area of the ASCA SIS ex- 
traction region (r = 4') is 6,400 times larger than the area 
of a Chandra extraction region (r = 3"). Moreover, the ra- 
dius of the ASCA SIS extraction region is twice as large as 
the radius of the entire optical galaxy. Thus any emission 
from hot plasma in the disk of NGC 4579 (analogous to 
the "Galactic Ridge" emission) is included in the ASCA 
spectrum and contributes to the strength of the Fe Ka 
line. The fact that the line energy was 6.7 keV rather 
than 6.4 keV bolsters this interpretation. To make a more 
quantitative assessment of the plausibility of our proposed 
explanation, we carried out the following tests: 

1. We extracted the total spectrum from the largest 
circular region that could reasonably fit in the Chan- 
dra field of view (centered on the AGN, with a ra- 
dius of 50 pixels or 24'.'5) and used it to look for an 
Fe Ka line, as described above. In this spectrum we 
detect the Fc Ka line at 68% confidence, with an 
equivalent width of 43 eV and we obtain an upper 
limit of 170 eV at 99% confidence, which is formally 
consistent with the ASCA measurement from 1998. 
This test illustrates that extended emission makes 
an increasing contribution to the equivalent width 
of the line in larger apertures. 

2. We used the spectral models for the AGN and the 
diffuse emission derived from the Chandra observa- 
tion to calculate the 2-10 keV spectrum and the 
apparent equivalent width of the Fe Ka line. For 
the purposes of this calculation we assumed the flux 
and spectrum of the AGN measured above and spec- 
trum of the extended emission determined in §5.2 
and summarized in Table 3, adopting Z = Zq and 
F x (2 - 10 keV) = 2.4 x 10" 13 erg s" 1 , which is 
apropriate for the entire circumnuclcar nebula. We 
find that the equivalent width of the Fe Ka complex 
in the combined AGN plus circumnuclear plasma 
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spectrum is approximately 40 eV. This implies that 
the entire strength of this line complex measured by 
ASCA can be accounted for, if we allow the entire 
galaxy to have 6 times the emission measure of the 
circumnuclear region, which is a reasonable possi- 
bility. 

6. VARIABILITY OF BRIGHT SOURCES 

To study the variability properties of the bright discrete 
sources in NGC 4736 and NGC 4579, we extracted their 
0.5-8.0 keV light curves from the same spatial regions used 
to extract spectra. Since only one quarter of the S3 CCD 
was read out during the observations of these two galax- 
ies, the frame time was 0.84 s. In the case of the AGN 
in NGC 4579 pileup is a cause for concern since it affects 
the properties of the light curve. The qualitative effect of 
pileup is to dilute the amplitude of any variability since it 
leads to some pairs of lower-energy photons being counted 
as single, higher-energy photons. Moreover, pileup also de- 
creases the S/N since some photon grades are transformed 
to unacceptable grades and those photons are not counted 
at all. The magnitude of these effects cannot be quantified 
without detailed and extensive simulations, which are well 
beyond the scope of this work. We have nevertheless made 
an effort to minimize the effects of pileup on the light curve 
of the AGN by restricting the band to 0.5-1.0 keV. This 
narrow, soft band encompasses the peak of the instrument 
sensitivity, thus including a significant fraction of the de- 
tected photons. The fact that it spans only a factor of 2 
in energy means that any pairs of lower-energy photons 
that are counted as one will not be included in the band. 
Contamination by pairs of photons below 0.5 keV counted 
within this band is small because the sensitivity of the 
CCD declines steeply below this energy. Hence, the most 
significant effect of pileup on light curves from this band 
is to reduce the S/N somewhat. Yet another way of min- 
imizing the effects of pileup on a light curve is to extract 
the light curve after excluding the core of the PSF, where 
pileup really occurs (cf, §5.4). Thus we have also extracted 
a light curve from an annulus with inner and outer radii 
of 0'.'5 and 4'.'2 respectively, which includes approximately 
half of the total counts. 

In Figure 7 we show the light curves of all the bright 
discrete sources in the two galaxies, binned into 30-minute 
time bins. These light curves show that the discrete 
sources vary on time scales of a few hours at the 20- 
30% level, while the AGN varies at the 5-10% level. The 
magnitude of the variability is quantified by applying the 
excess variance test (Nandra et al. 1997a) to these light 
curves, whose results are summarized in Table 5. In the 
case of NGC 4579, the excess variance computed from the 
soft light curve agrees with the value obtained from the 
outer-PSF light curve and is a factor of 2 higher than the 
full-band light curve. This suggests that pileup dilutes 
the variability in the full-band light curve but not in the 
other two light curves. It also suggests that the shape 
of the spectrum does not change as the light curve flick- 
ers. The level of variability of the AGN in NGC 4579 is 
higher than what is observed in typical LINERs of com- 
parable luminosity (Ptak et al. 1998) and approaches that 
of Seyfert galaxies whose luminosities are 1-2 orders of 
magnitude higher than that of NGC 4579. Unfortunately, 



it is difficult to draw more extensive conclusions from the 
variability properties of the nucleus of NGC 4579 because 
of the caveats associated with pileup. 

We have also searched for coherent, periodic signals in 
the light curves of the bright discrete sources, but no such 
signals were found. We carried out this search by com- 
puting the Fourier transforms of the source light curves 
up to frequencies of 0.6 Hz, following the methodology of 
Eraclcous, Patterson, & Halpern (1991). Using the noise 
properties of the power spectra we have estimated upper 
limits to the amplitudes of undetected signals, which we 
list in Table 5. These limits are expressed as a minimum 
detectable "pulsed fraction," which is the amplitude of a 
sinusoidal signal relative to the mean level of the light 
curve. 

7. DETECTION AND POPULATION PROPERTIES OF THE 
DISCRETE X-RAY SOURCES 

To study the population of discrete X-ray sources in 
the target galaxies, we used the source detection program 
WAVDETECT (Freeman et al. 2001), which is included in the 
CIAD package. We searched for sources on several spatial 
scales (from 1.4 to 11.6 pixels, in steps of y/2), accepting 
detections associated with less than 10~ 6 chance probabil- 
ities due to local background fluctuations. We repeated 
the search twice, first in the full, 0.5-8 keV, band and 
then in the soft, 0.5-2 keV, band. We combined the re- 
sults of the two passes to compute an X-ray color for each 
source, which we define as the ratio of counts in the hard 
(2-8 keV) band to counts in the soft (0.5-2 keV) band. 
Because the sensitivity of the instrument rises steeply to- 
wards low energies, all of the sources are detected in the 
total and soft bands, but nearly half of the sources have a 
hard band count rate consistent with zero. In the end our 
detection limit corresponds to about 6 counts per source, 
although it varies over the field of view because of (a) vi- 
gnetting, (b) degradation of the image quality off axis, and 
(c) patchy diffuse emission, which increases the effective lo- 
cal background. After carrying out the above automatic 
procedure, we checked the measured properties of many 
of the detected sources manually. We paid particular at- 
tention to sources in the crowded field in the nucleus of 
NGC 4736 and we excluded any sources falling along the 
readout streak in the image NGC 4579. We tabulate the 
properties of detected sources in Tables 6 and 7, where we 
give the coordinates of each source, the total count rate, 
the X-ray color, and an estimate of the 0.5-8 keV luminos- 
ity. For sources whose hard count rate is consistent with 
zero, we only give an upper limit to the color. 

To estimate the luminosities given in the source tables 
we assumed that the source spectra could be described by 
a simple power law with an absorbing column equal to 
the Galactic column. With this assumption we simulated 
spectra with a wide range of photon indices (0.5-5.0) and 
measured and tabulated their colors and count rates in the 
same way as we did for the observed spectra. We used the 
resulting lookup table to estimate the photon indices and 
fluxes of detected sources from their observed colors and 
count rates. By comparing the results of this method with 
results obtained for bright sources for which spectra were 
extracted and fitted with models, we estimate the uncer- 
tainty in the photon index and flux to be about 10-20%. 
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This uncertainty is in addition to the Poisson uncertainty 
on the count rate listed in Tables 6 and 7. In cases when 
only an upper limit to the color was determined from the 
data, we used this upper limit to determine the photon in- 
dex. Thus the uncertainty in the flux of these sources can 
be up to 50%, judging from the relation between photon 
indices and count rates in simulated spectra. 

As we already mentioned in §4, in NGC 404 we do not 
find any X-ray sources other than the nucleus, down to a 
limiting luminosity of 1 x 10 36 erg s _1 . In NGC 4736 and 
NGC 4579 we find respectively 39 and 21 sources, with 
luminosities in the range 1 x 10 36 — 1 x 10 39 erg s _1 in the 
former galaxy and 1 x 10 37 — 9 x 10 39 erg s _1 in the lat- 
ter. The lower end of the range is effectively our detection 
limit. Thus, in the case of NGC 4736 we can detect typi- 
cal supernova remnants (SNRs) and Be-star X-ray binaries 
(BcXBs), while in the case of NGC 4579, we can detect 
low-mass X-ray binaries (LMXBs) or more luminous ob- 
jects. We have estimated the probability of contamination 
of the source population by unrelated, background sources 
using the number counts of sources derived from deep 
Chandra surveys by Tozzi et al (2001). In the observation 
of NGC 4736 the limiting flux is 4.8 x 10" 16 erg cur 2 s _1 , 
which yields an expectation value of 5 background sources 
within the field of view, while in the observation of 
NGC 4579 the limiting flux is 5.6 x 10~ 16 erg cm" 2 s" 1 
and the expected background source number is 4. These 
are conservative upper limits to the number of background 
sources, since the flux of such sources would be attenuated 
as it passes through the disk of the foreground galaxy. If 
we assume, for example, that the column density perpen- 
dicular to the galactic disks is 1 x 10 21 cm~ 2 , the expected 
number of background sources in each case decreases by 1 . 
Thus, contamination from background sources is negligible 
from the point of view of our statistical studies, although 
it may be an issue in our discussion of the properties of 
individual sources. We also considered the possibility of 
foreground contamination, i.e. that some of the detected 
X-ray sources are foreground stars. The absence of op- 
tical counterparts of any of the detected X-ray sources 
(see §4) indicates that there is no foreground contamina- 
tion. This is not unexpected since the Galactic lattitude of 
the target galaxies is rather high (—27°, +76°, and +74° 
for NGC 404, NGC 4736, and NGC 4579 respectively). 
The faintest optical objects that could be associated with 
detected X-ray sources are M stars. Using the X-ray-to- 
optical flux ratios of M stars reported by Hempellmann 
et al. (1995), we estimate that M star counterparts to the 
sources that we have detected should have R < 20.2, and 
would have been included in the USNO A2.0 catalog. 

The spatial distribution of sources differs markedly be- 
tween NGC 4736 and NGC 4579. In the former galaxy 
there is a concentration of sources within 10" from the 
nucleus, and a separate concentration at the radius of the 
star-forming ring (~ 45-50"). Thus, the distribution of 
X-ray sources does not follow the distribution of infrared 
light (i.e., starlight) from the host galaxy, but rather the 
distribution of emission-line light (cf, Smith et al. 1994), 
suggesting a connection with star formation. In contrast, 
the sources in NGC 4579 are distributed smoothly with 
distance from the center. 

In Figure 8 we show the "color-magnitude" diagrams 



and the luminosity distributions of the source populations 
of the two galaxies. In both galaxies, the most luminous 
sources have photon indices between 1.0 and 2.0. At lower 
luminosities only limits on the photon index are avail- 
able, although there are some notable exceptions in both 
galaxies, where we detect rather hard sources (r ~ 1) at 
Lx ~ 10 37 erg s" 1 . The luminosity distribution appears 
to peak just below 10 37 erg s^ 1 in NGC 4736 and just 
below 10 38 erg s" 1 in NGC 4579. 

In Figure 9 we compare the cumulative X-ray source 
luminosity functions of NGC 4736 and NGC 4579 with 
those of other galaxies with published Chandra data. We 
include two early- type galaxies, NGC 4697 and NGC 1553 
(E and SO, respectively; Sarazin et al. 2001; Blanton et al. 
2001) and two starbursting spirals, M82 and the Circinus 
galaxy (Griffiths et al. 2000; Bauer et al. 2001). NGC 4579, 
NGC 4697, and NGC 1553 stand out in this figure be- 
cause their cumulative X-ray source luminosity functions 
are steeper than those of the other three galaxies. This 
difference can be quantified by fitting a simple power-law 
model [N(L) cx L~ a or N(> L) cx L"^" 1 )] to the data 
at luminosities above the initial plateau. In the case of 
NGC 4579, NGC 4697, and NGC 1553 we find power- 
law indices (a) between 2.0 and 2.4, while in the other 
two galaxies the power-law indices range between 1.4 and 
1.6. In fact, the difference is more extreme than what the 
above power-law indices indicate because the X-ray lumi- 
nosity functions of the two early-type galaxies steepen at 
luminosities higher than about 10 38 erg s -1 (Sarazin et al. 
2001; Blanton et al. 2001). 

It is interesting to note that the galaxies with steep X- 
ray source luminosity functions are the ones without signs 
of current or recent star formation activity. In this con- 
text it is noteworthy that a similar effect is observed in two 
nearby spiral galaxies, M31 and M81 (Primini, Forman, & 
Jones 1993; Tennant et al. 2001), where the X-ray luminos- 
ity functions of bulge and disk sources have been measured 
separately. The bulge luminosity functions resemble those 
of early type galaxies, while the luminosity function of the 
disk of M81 resembles that of the starbursting spirals (the 
disk luminosity function of M31 has not been reported). 
In §8 we interpret this difference in X-ray source luminos- 
ity functions by appealing to current or recent intense star 
formation and the age of the underlying stellar population. 

8. DISCUSSION 

8.1. Implications for the Central Engines of LINERs 

Our original aim in undertaking this program was to 
search for the power source of the optical emission lines, 
which are the hallmark of LINERs. Our main result is 
that the objects in our small collection arc quite diverse 
in their properties, with a different power source appar- 
ently playing the dominant role in each galaxy. As a first 
step in considering the implications of these findings, we 
compare them with published results of previous X-ray 
observations of these galaxies with ASCA. In Table 8 we 
list the 0.5-10 keV luminosity of the different components 
contributing to the X-ray emission from each galaxy. This 
table quantifies the general description given in §4. The 
sum of these components reproduces fairly closely the flux 
detected by ASCA from each of these galaxies. One no- 
table discrepancy between the Chandra and ASCA spectra 
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X-ray luminosity and less than half of the UV luminos- 
ity. 6 The properties of the non-thermal radio continuum 
emission from the nucleus (Turner & Ho 1994) are con- 
sistent with an origin in a group of several, young SNRs 
of the Cas A variety, whose kinetic energy heats the gas 
responsible for the diffuse X-ray emission (cf, the discus- 
sion of the diffuse X-ray emission of the LMC by Wang 
et al. 1991). In the same spirit, the nuclear, point-like, 
UV source is very likely an OB association, which hosts 
a HMXB (NGC 4736 X-2). The observations of Maoz et 
al. (1995) provide a lower limit to the monochromatic lu- 
minosity at 2270 A (the central point source is saturated) 
of i A (2270 A) > 4.4 x 10 36 erg s" 1 A" 1 , which is equiva- 
lent to the luminosities of 5 late-0 supergiants. Indepen- 
dent evidence from optical spectroscopy by Taniguchi et 
al. (1996) leads to similar conclusions about the nature of 
the activity in the nulceus of NGC 4736. These authors 
find that the optical absorption-line spectrum of the nu- 
cleus of NGC 4736 agrees with models of a 10 9 -year old 
starburst. They also find that the nuclear Ha luminosity 
is no larger than that of luminous Galactic H II regions 
and that six 06 stars are enough to power it. Thus, this 
particular LINER seems to be associated with a recent or 
aging starburst rather than an AGN. 

8.2. The Discrete Source Populations 

The discrete sources detected in NGC 4736 and 
NGC 4579 may include a variety of objects. The most 
luminous sources (Lx > 10 37 erg s _1 ) are most likely X- 
ray binaries: a combination of LMXBs and HMXBs. This 
conclusion is supported by two pieces of evidence: their 
spectra, which are similar to those of bright Galactic X- 
ray binaries, and their variability properties, which suggest 
that these are individual objects, rather than clusters of 
less luminous sources. Other types of object, such as SNRs 
and BeXBs, could contribute to the population of lower lu- 
minosity sources, detected primarily in NGC 4736. If we 
could determine the colors of the lower-luminosity sources 
we could perhaps use them to differentiate between SNRs 
and BeXBs since the 2-10 keV spectra of the former are 
dominated by line emission at low energies, while the spec- 
tra of the latter resemble power laws. 

The high concentration of luminous sources in the nu- 
cleus of NGC 4736 suggests that we are observing the af- 
termath of an intense episode of star formation. By ex- 
tension, a significant fraction of the discrete sources in 
NGC 4736 could be HMXBs rather than LMXBs. Mod- 
els of populations of interacting binaries associated with 
bursts of star formation (Cerviho, Mas-Hesse, & Kunth 
1997; Van Bever & Vanbeveren 2000; Sipior & Eracleous 
2000) show an initial hard-X-ray luminous phase, which is 
intimately connected with the rise and fall of HMXBs (see 
also the discussion of the processes by Dalton & Sarazin 
1995). The first HMXBs make their appearance after the 
death of the most massive O stars, and shine brightly until 
the demise of their O and B secondaries. At this point the 
total hard X-ray luminosity of the population should de- 
cline until the first luminous LMXBs arrive on the scene, 
causing the X-ray luminosity to level off with time. Thus, 
the hard X-ray emission from the population should per- 

6 For the sake of completeness, we note that the stellar velocity dispersion of NGC 4736 of 136 km s _1 (Mc Elroy 1995) implies a black hole 
mass of 2 — 3 X 10 7 Mq and an Eddington ratio of L/L^dd ~ 3 X 10~ 6 (see §8.4 for the details of the methodology). 



is the detection of Fe Ken lines at 6.4 and 6.7 keV in the 
ASCA spectra of both NGC 4736 and NGC 4579 (Roberts 
et al. 1999; Terashima et al. 1998). Such lines arc not seen 
in the Chandra spectra of any of the bright discrete sources 
in these galaxies. What Chandra does detect, however, is 
a hot component (kT ~ several keV) to the plasma that 
permeates the central regions of these two galaxies. The 
discrepancy can thus be resolved if the Fe Ka emission 
lines in the ASCA spectra arise from this hot plasma, as 
discussed in §5.4. 

In NGC 4579 the case for an AGN seems iron clad. A 
central, compact source dominates the power budget, as 
we illustrate in Table 8. The spectral energy distribution 
of this source has a significant contribution from UV light, 
broad optical and UV emission lines, as well as a non- 
thermal radio spectrum with a brightness temperature of 
2 x 10 8 K (Falcke et al. 2000). All these properties are 
classic traits of an AGN. 

At the opposite extreme, NGC 404 harbors a rather 
anemic star- forming region in its nucleus. Its X-ray lu- 
minosity of 1 x 10 37 erg s _1 is comparable to that of a 
single HMXB or a giant star-forming region such as 30 
Doradus (cf, Wang & Helfand 1991). The soft, thermal 
X-ray spectrum favors an origin of the X-rays in the hot 
gas of a superbubble, powered primarily by off-center su- 
pernovae (see the discussion by Chu & Mac Low 1990). 
This view is supported by the absorption-line UV spec- 
trum and the monochromatic luminosity at 2270 A of 
L A (2270 A) = 1.2 x 10 36 erg s" 1 A" 1 (from Maoz et 
al. 1995), which is equivalent to that of 10 late-0 su- 
pergiants. The sound speed in the hot gas is of order 
c s w 380 (fcT/0.5 keV) 1/2 km s _1 , which is considerably 
larger than the velocity dispersion of the stars in the nu- 
cleus of the galaxy (<r* = 55 km s _1 ; Faber et al. 1989). 
Thus the gas must be unbound, a conclusion bolstered by 
the "blow-out" morphology of the Ha+[N II] image, which 
also lends support to the superbubble hypothesis for the 
X-ray emitting gas. 

In the case of NGC 4736, there is no compelling evi- 
dence for the presence of an AGN even though this galaxy 
is fairly bright in X-rays and harbors an unresolved UV 
source. The X-rays are produced by a dense cluster of 
stellar X-ray sources in the inner 2 kpc of the galaxy, and 
to a lesser degree by diffuse emission from hot gas (see the 
energy budget in Table 8). This possibility was correctly 
anticipated and discussed by Roberts et al. (1999) and 
Roberts et al. (2001) who nonetheless preferred an AGN 
interpretation in the end. As we argue below, these X- 
ray sources could have been produced in a recent episode 
of star formation and probably consist of high-mass X- 
ray binaries (HMXBs), LMXBs,and SNRs. The bright- 
est nuclear X-ray sources have very similar X-ray spec- 
tra and comparable luminosities, making the association 
of any one of them with an AGN rather implausible. It 
should be pointed out, however, that source X-2, coin- 
cides with the nuclear UV source detected by the HST; 
thus it stands out among the nuclear X-ray sources as 
the only one with a UV counterpart. But even if we as- 
sociate this source with an AGN, it is not very impor- 
tant energetically since it contributes less than 20% of the 



10 



Eracleous et al. 



sist for up to 50 Myr after the vigorous star formation has 
stopped. Inspired by the model predictions, we propose 
that the nucleus of NGC 4736 is an "aging" starburst or 
in an immediate post-starburst phase, where the HMXBs 
dominate the top of the X-ray source luminosity function. 
In support of our hypothesis we note the qualitative dif- 
ference between the cumulative luminosity functions of X- 
ray sources in nearby galaxies, shown in Figure 9 and dis- 
cussed in §7. Galaxies currently undergoing vigorous star 
formation, such as M82 and Circinus, have cumulative lu- 
minosity functions that are flatter than those of galaxies 
with primarily old stellar populations such as NGC 4697 
and NGC 1553. We speculate that the flatter slopes re- 
sult from an excess of luminous HMXBs that have formed 
recently. The steep luminosity function of NGC 4579 fits 
well within this picture since there are no indications of 
ongoing, intensive star formation in its bulge or disk. 

A particularly striking source is NGC 4579 X-l, which 
could be a UXB. If this object is at the distance of 
NGC 4579, as we have been implicitly assuming, its 2- 
10 keV luminosity would exceed 10 40 erg s _1 . Its spec- 
trum is very similar to the spectra of the X-ray binaries 
in NGC 4736 and, in fact, it is the only one of the bright 
X-ray sources whose spectrum is better described by a 
Comptonized blackbody model than by a simple power 
law. This object is consequently very similar to luminous 
Galactic X-ray binaries. Its luminosity corresponds to the 
Eddington limit for a 70 M© object, which could be re- 
garded as a lower limit on the mass of the accreting com- 
pact object, at first glance. However, the Eddington limit 
need not apply in the non-spherical accretion geometry of 
X-ray binaries, and a number of cases are known in which 
it is blatantly violated. There are at least 3 X-ray pul- 
sars in the Magellanic Clouds (LMC X-4, SMC X-l, and 
A 0538-66; Nagase 1989; Woo, Clark, & Levine 1995, and 
references therein) whose luminosities exceed the Edding- 
ton limit for a neutron star by as much as a factor of 
5 (and even by a factor of 30 in the flares of LMC X-4; 
Woo ct al. 1995). Of course violations of the Eddington 
limit may be transient (cf, the discussion of SMC X-l by 
Hclfand & Moran 2001) or merely appear as such because 
of beaming (King et al. 2001). However, the fact remains 
that at certain times or from certain viewing directions 
some X-ray binaries appear extremely luminous. Thus, 
the accreting compact object in NGC 4579 X-l need not 
be more massive than about 10 M©. Another possibility 
is that this source is at a large distance behind NGC 4579. 
The fact that we expect 3 background sources in the image 
of NGC 4579, as well as the large measured column den- 
sity (comparable to the column density perpendicular to 
the disk of a spiral galaxy) make this a plausible scenario. 
It should be pointed out, however, that the 3 expected 
background sources should be rather faint according to 
the flux distributions measured in deep surveys. The a pri- 
ori probability of finding a background object with the 2— 
10 keV flux of NGC 4579 X-l in our field of view is approx- 
imately 6%. If this source were a background AGN, then 
its redshift, as inferred from the observed flux, would be 

z = 0.01 (ffo/50 km s" 1 Mpc" 1 ) (i x /10 42 erg s" 1 ) 172 , 
suggesting that the source could easily be a garden variety 
Seyfert galaxy at z = 0.01 or a luminous quasar at z = 1. 
Ascertaining the nature of this source will have to await 



identification of its optical counterpart. 

8.3. The AGN and Extended Emission in NGC 4579 

Using the uncontaminatcd luminosity of the AGN in 
NGC 4579, we can explore the properties of the accretion 
flow onto the supcrmassive black hole that powers it. We 
can estimate the mass of the central black hole using re- 
cently established correlations between it and the velocity 
dispersion of stars in its immediate vicinity (Gebhardt et 
al. 2000; Merrit & Ferrarese 2001). The velocity disper- 
sion of 185 km s _1 measured by Dressier (1987) yields a 
black hole mass in the range 6 to 13 x 10 M Q , compa- 
rable to what is found in Seyfert galaxies (Ferrarese et 
al. 2001). Assuming that the observed 0.5-10 keV lu- 
minosity is about 10% of the bolometric luminosity, we 
can estimate the Eddington ratio for this AGN to be 
L/LEdd ~ 1 — 2 x 10~ 4 , which is a low value, consistent 
with expectations for a low-luminosity AGN. At such a 
low Eddington ratio, the inner accretion flow is expected 
to be advection-dominated (an ADAF or similar structure) 
and thus very hot, vertically extended, and optically thin. 
Such a structure would be very different from a conven- 
tional, optically thick, geometrically thin accretion disk 
in its observational signature. In particular, the absence 
of an optically thick inner accretion disk can explain the 
weakness of the Fe Ka line of NGC 4579 relative to what 
is observed in Seyfert galaxies (see the detailed discussion 
in Eracleous, Sambruna, & Mushotzky 2000). A vertically 
extended ADAF comprising the inner accretion structure 
can illuminate the outer thin disk and drive the emission 
of optical Balmer lines as well as X-ray Fe Ka lines. This 
scenario offers an appealing explanation for the double- 
shouldered profile of the Ha line of NGC 4579 (Barth ct 
al. 2001) just as it does for the double-peaked Balmer 
lines in broad-line radio galaxies (Eracleous & Halpern 
1994). However, models for the continuum spectra of 
ADAFs (e.g., Narayan, Yi, & Mahadevan 1996) predict 
X-ray spectral indices flatter than the typical values ob- 
served in Seyfert galaxies. This does not seem to hold for 
NGC 4579, where the 2-10 keV spectral index is typical 
of a Seyfert galaxy but the equivalent width of the Fe Ka 
line is not. However, the UV-X-ray spectral index (a ox ) 
of NGC 4579 is rather typical of LINERs as shown by Ho 
(1999). Interestingly enough, the variability amplitude is 
also high for a low-luminosity AGN and is comparable to 
what is observed in Seyferts (see §6). Therefore, NGC 4579 
appears to be a "hybrid" object with some properties char- 
acterisitc of Seyferts and other characteristic of LINERs 
or low-luminosity AGNs. 

The heating mechanism of the extended emission region 
of NGC 4579 is unclear. We can confidently exclude any 
mechanical interaction between it and an AGN jet since 
the VLBI radio images of Falcke et al. (2000) show only 
an unresolved radio core with no extension. One possible 
mechanism is that the extended emission region is pho- 
toionized by the AGN. The X-ray luminosity of the AGN 
is 20 times higher than that of the circumnuclear neb- 
ula and a fraction of the hard X-ray luminosity from the 
circumnuclear region may be contributed by BeXBs and 
LMXBs, which we have not detected. Thus, photoiniza- 
tion may be a viable mechanism. Moreover, as Maoz et 
al. (1995) point out, the extrapolation of the observed UV 
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luminosity is enough to power the optical emission lines 
from the same region. Another mechanism is heating by 
stellar processes in the wake of a circumnuclear starburst. 
In the context of this scenario, the starburst should be old 
enough that its luminous O stars have passed away and 
supernova remnants have blended together into a diffuse 
circumnuclear medium. Unfortunately, we cannot assess 
the likelihood of the above scenarios with the information 
currently available. 

9. EPILOGUE 

The results of our exploratory observations indicate that 
LINERs are a heterogeneous population of objects. This 
should not come as a surprise, in view of the way that the 
class is defined. But an important question that remains 
outstanding is what fraction of LINERs are powered by 
accretion onto a supermassive black hole and how these 
objects can be identified unambiguously. Unfortunately, 
the first part of this question cannot be answered until the 
second part is answered first. One rather clear signature 
of an AGN is the presence of broad emission lines in the 
optical or UV spectra of LINERs. However, the absence 
of broad emission lines cannot be taken as evidence that 
there is no AGN. This is especially true when one considers 
the practical difficulties of finding weak broad lines in the 
starlight-dominated spectra of LINERs. Another possibly 
clear signature of an AGN is the presence of a non-thermal 
radio core with a high brightness temperature. 

X-ray observations can also shed light on this question, 
although they are more demanding in terms of resources 
and may not always yield a definitive answer. The high 



spatial resolution of Chandra now allows us to examine 
the morphology of LINER nuclei in detail and combine 
morphological and spectral information, just as we have 
done here. However, one can also envisage cases where a 
LINER harbors a single, discrete source of "intermediate" 
luminosity (10 39 -10 40 erg s _1 ), which could equally well 
be an AGN or a UXB. Moreover, observations with ASCA 
should be regarded with caution, since they can sometimes 
be misleading because of their poor spatial resolution. A 
good case in point is provided by NGC 4736, whose X-ray 
binaries and diffuse emission produce a spatially integrated 
spectrum that can be modeled by a hard power law and 
a thermal plasma component. This "canonical" LINER 
spectrum can result from a current or recent starburst, 
just as well as it can result from an AGN. 
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Table 1 

Targets and Their Basic Properties 









Spatial 


Galactic 


Galactic 




Target 


Hubble 


Distance b 


Scale 


A H C 


E(B — V) d 




Name(s) 


Type a 


(Mpc) 


(PC/") 


(10 20 cm- 2 ) 


(mag) 


(km s x ) 


NGC 404 


SA(s)0-: 


2.4 


12 


5.30 


0.059 


-46 


NGC 4736 (M94) 


(R)SA(r)ab 


4.3 


21 


1.41 


0.018 


307 


NGC 4579 (M58) 


SAB(rs)b 


16.8 


81 


2.54 


0.041 


1521 



References. - (a) de Vaucouleurs et al. (1991), (b) Tully (1988), (c) Stark et al. (1992), (d) 
Schlegel, Finkbeiner, & Davis (1998), (c) Ho et al. (1997a) 



Table 2 
Observation Log 





Observation 


Exposure 




Start 


Time 


Target 


(UT) 


(s) 


NGC 404 


1999 Dec 19, 21:12:43 


23,864 


NGC 4736 


2000 May 13, 03:12:11 


47,366 


NGC 4579 


2000 May 02, 16:26:44 


33,905 



Table 3 

Parameters of Thermal Plasma Models 











NGC 4579 6 


Model Parameters 


NGC 404 a 


NGC 4736 b 


Z =free 


Z = Z Q 


A H (10 20 cm" 2 ) 


< 6 


1.4 C 


< 6 


7+9 
'-4 


Abundance (Z/Zq) 


< 0.12 


3+ 1 - 4 
u.o_ 2 


2+ 1 ' 3 


1 (fixed) 


kT d (keV) 


fi +0 - 3 

U.D_ 4 


0.56 ±0.06 


0.3 ±0.2, 0.6±° i, 


> fi ii+ - 05 n ^8 +0 - 04 10+ 44 

> U.ll_ 04 , U.OO_ 08 , 1U_4 


EM e (10 61 cm" 3 ) 
Total x 2 


0.35 


0.7 


33, 34, 28 


34, 11, 27 


10.45 


3.159 


22.507 


25.546 


Reduced x 2 


0.696 


0.632 


1.072 


1.161 


Flux (erg cm -2 s _1 ) 










0.5-2 keV 


1.5 x 10- 14 


2.0 x 10- 14 


1.2 x 10- 13 




2-10 keV 






1.2 x 10- 13 




Luminosity (erg s- 1 ) 










0.5-2 keV 


1.0 x 10 37 


4.4 x 10 37 


4.1 x 10 39 




2-10 keV 






4.1 x 10 39 





a The spectrum was extracted from a circular region of area 13 square arcseconds, equal to the area of a 
point source. 

b The spectrum was extracted from a region of area 239 square arcseconds. 
c Held fixed at the Galactic value because it was unconstrained. 

d Single-temperature models for NGC 404 and NGC 4736, 3-temperaturc model for NGC 4579. 
e Emission measure: EM = J n c rii dV, where n c and ri; are the electron and ion densities. 
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Table 4 

Power-Law Fits to Spectra of Bright Discrete Sources 





Column 


Power-Law 


Observed 


2-10 keV 








Density 


Photon 


2-10 keV Flux 


Luminosity a 


Total 


Reduced 


Object 


(10 20 cm- 2 ) 


Index 


(erg cm -2 s -1 ) 


(erg s" 1 ) 


x 2 


x 2 


NGC 4736 














X-l 


< 3.7 


1.13 ±0.10 


6.4 x 10~ 13 
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34.13 


0.7937 
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30.19 


1.078 
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1.6 ±0.2 
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0.9234 


NGC 4579 














X-l 
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1.8 ±0.2 


2.7 x 10~ 13 


9.1 x 10 39 


22.05 


0.9588 



a Intrinsic luminosity corrected for absorption. 



Table 5 

Variability Properties of Bright X-Ray Sources 



Object 


Excess Variance Maximum Pulsed Fraction 

(s- 1 ) (%) 


NGC 4736 






X-l 


0.06 ±0.02 


13 


X-2 


0.06 ±0.04 


10 


X-3 


0.10 ±0.07 


17 


X-4 


< 0.1 


19 


NGC 4579 






AGN (0.5-1 keV) 


0.035 ± 0.008 




AGN (0.5-8 keV) 


0.017 ±0.001 




AGN (0.5-8 keV), PSF wings 


0.05 ±0.02 




X-l (UXB) 


< 0.02 
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Table 6 

Catalog of Discrete Sources in NGC 4736 



Official Name 
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Table 7 

Catalog of Discrete Sources in NGC 4579 



Official Name 
(CXOU J) 


Alias 




Coordinates 
(J2000) 




Count Rate 
(10~ 2 s- 1 ) 


Color 
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Table 8 
Power Budget 



0.5-10 keV Luminosity (erg s x ) 



NGC 404 NGC 4736 NGC 4579 

AGN 

Diffuse or Thermal Emission a 
UXBs 

Other discrete sources 



Chandra Total 

ASCA Total b < 5 x 10 37 3 x 10 39 3.0 x 10 



a This is the total luminosity of the entire diffuse source, estimated after 
subtracting the contribution of point sources embedded in it. 

b References. - NGC 404: Roberts et al. (2001); NGC 4736: Roberts et 
al. (1999); NGC 4579: Terashima et al. (1998) 



2.9 x 10 41 

1 x 10 37 6 x 10 38 1.6 x 10 40 
... 1.3 x 10 40 
3.7 x 10 39 2 x 10 39 



1 x 10 37 4 x 10 39 3.2 x 10 
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Fig. 2. — Montage of Chandra images of the target galaxies. For each galaxy we show an image encompassing all the interesting features in 
the field of view as well as an enlarged view of the nucleus. All but one of the frames are labeled with the right ascension along the horizontal 
axis and the declination along the vertical axis (the coordinate epoch is J2000, north is up, and east is to the left). The image of the nucleus 
of NGC 404 has been deconvolved using the Lucy-Richardson algorithm; its frame is calibrated in pixels, where each pixel corresponds to 
0'.'15. The discrete "blobs" around the nucleus of NGC 404 are not discernible in the original, unprocessed image. A horizontal bar in each 
image gives the linear scale at the distance of that galaxy. Bright X-ray sources, referred to in the text and in the tables, are labeled on the 
images. 
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1 2 3 0.5 1 1.5 3 2 4 

Energy (keV) Energy (keV) Energy (keV) 



Fig. 3. — Spectra of the thermal plasma regions in the three target galaxies, the nucleus of NGC 404 and the extended, circumnuclcar 
emission in NGC 4736 and NGC 4579. For each object we show the observed spectrum with the best-fitting model superposed (upper panel) 
and the post-fit residuals scaled by the error bar (lower panel). The models consist of linear combinations of thermal plasma components, with 
abundances left as free parameters (1 component in NGC 404 and NGC 4736, and 3 components in NGC 4579). The parameters describing 
these models arc summarized in Table 3. 




Fig. 4. — Spectra of the brightest X-ray binaries in NGC 4736 and NGC 4579. The objects are labeled according to the conventions defined 
in Tables 6 and 7. For each object we show the observed spectrum with the best-fitting model superposed (upper panel) and the post-fit 
residuals scaled by the error bar (lower panel). The parameters describing the best-fitting models are summarized in Table 4. 
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Energy (keV) 



Fig. 5. — Spectrum of the AGN in NGC 4579 with the best-fitting model superposed. The model consists of a simple power law modified 
by photoelectric absorption in the Galaxy, as described in §5.4. The lower panel shows the residuals scaled by the error bar at each bin. The 
large residuals around 0.6 keV may be the result of poor calibration of the S3 CCD. Due to calibration uncertainties, we are not able to assess 
the reality of this feature. 
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Fig. 6. — Confidence contours showing the 99% upper limits on the equivalent width of the Fe Ka line in the spectrum of the AGN in 
NGC 4579 as a function of its FWHM. The equivalent width scale has been corrected for the effect of pileup by multiplying all values by 
1 + /pi where / p is the pileup fraction in the 6-7 keV band (see §5.4 of the test for a detailed discussion). The three curves refer to different 
assumed line energies; solid curve: 6.4 keV, dashed curve: 6.7 keV, dot-dashed curve: 6.9 keV. The error boxes (dotted lines corresponding 
to the 90% confidence regions) show the measurements from ASCA observations (Terashima et al. 1998, 2000c) with dates as marked. 
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Fig. 7. — Light curves of bright point sources in NGC 4736 and NGC 4579. The objects are labeled according to the conventions defined 
in Tables 6 and 7. In the case of the AGN in NGC 4579, we plot the 0.5-1.0 keV and 0.5-8.0 keV light curves, as well as the 0.5-8.0 keV 
light curve extracted from the wings of the PSF. In all other cases we plot only the 0.5-8.0 keV light curves. 
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NGC 4736 Color-Magnitude Diagram 
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NGC 4579 Color-Magnitude Diagram 
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Fig. 8. — Color-Magnitude diagrams and luminosity distributions for the X-Ray binaries in NGC 4579 and NGC 4736. In the top frame of 
each set the vertical axis is calibrated in color (left) and in spectral index (right). Upper limits to the color are denoted by arrows. 
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Fig. 9. — Cumulative X-ray source luminosity functions of NGC 4579 (dot-dashed line) and NGC 4736 (solid line) compared to those of 
other galaxies observed by Chandra. The comparison galaxies arc NGC 4697 and NGC 1553 (E and SO, dotted lines), and M82 and Circinus 
(starbursting spirals, solid lines). The UXB in NGC 4579 is excluded from its luminosity function in the interest of clarity. 



